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Weak gravitational lensing by intervening large-scale structure induces a distinct signature in the 
cosmic microwave background (CMB) that can be used to reconstruct the weak-lensing displacement 
map. Estimators for individual Fourier modes of this map can be combined to produce an estimator 
for the lensing-potenial power spectrum. The naive estimator for this quantity will be biased upwards 
by the uncertainty associated with reconstructing individual modes; we present an iterative scheme 
for removing this bias. The variance and covariance of the lensing-potenial power spectrum estimator 
are calculated and evaluated numerically in a ACDM universe for Planck and future polarization- 
sensitive CMB experiments. 
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_ ■ The primordial cosmic microwave background (CMB) was generated when photons first decoupled from the baryonic 
fluid when the universe was only 400,000 years old. The vast majority of these photons travel unperturbed to the 
present day, and features of their angular power spectrum such as acoustic peaks and the damping tail [1] record 
valuable information about cosmological parameters [2] . Baryons and dark matter evolve from small inhomogeneities 
at decoupling into increasingly complicated large-scale structure which can subtly perturb the observed pattern of 
CMB anisotropies. Assuming that the primordial CMB is Gaussian, non-Gaussian correlations in the observed map 

■ can be used to reconstruct the intervening large-scale structure [3] . In addition to the importance of learning about 

■ the large-scale structure itself, reconstruction of the weak-lensing potential generated by structure is essential to 
, constraining tensor perturbations. Weak lensing converts a fraction of the E-mode polarization generated by scalar 

perturbations at the last-scattering surface into B-mode polarization in the observed map. Only by subtracting this 
B-mode polarization can one conclusively detect the primordial B-modes which serve as a model-independent signal 
of tensor perturbations [4]. Understanding lensing reconstruction requires a more detailed discussion of how weak 
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' lensing affects the CMB. 

Weak gravitational lensing deflects the paths of CMB photons as they travel from the last-scattering surface to the 
observer. This deflection is accomplished by a projected lensing potential which is a weighted line-of-sight integral of 
the gravitational potential between the observer and the surface of last-scattering. At each point on the sky, lensing 
remaps the temperature and polarization to that of a nearby point at the last-scattering surface, the deflection angle 
being the gradient of the aforementioned projected lensing potential. Assuming that this deflection angle is small, 
the temperature at any point can be expanded in a Taylor series in the gradient of the lensing potential. In Fourier 
space, this expansion appears as a series of convolutions of individual temperature and projected potential modes. 
The observed temperature-squared map in Fourier space also appears as a convolution of individual Fourier modes. 
Subject to an overall normalization dependent on the scale of the Fourier mode, these convolutions cancel in such a 
manner that each Fourier mode of the temperature-squared map acts as an estimator for the same Fourier mode of 
jj] ■ the projected lensing potential. 

Lensing reconstruction as outlined above has been considered previously [5,6]. In these works, two sources of noise 
were identified, and a filter of the temperature-squared map in Fourier space was chosen to minimize the variance 
associated with lensing reconstruction subject to these noise sources. The first source is intrinsic signal variance; the 
observed large-scale structure is one arbitrary member of an ensemble of realizations allowed by theory. The second 
source of noise, endemic to this method of lensing reconstruction, is a consequence of the nature of the primordial 
CMB. Like the large-scale structure itself, the pattern of CMB anisotropies at the last-scattering surface is only one 
of many possible realizations allowed by theory. We do not know a priori which of these realizations nature has 
provided us, and this uncertainty hinders our ability to deconvolve the effects of lensing from true anisotropies at 
the last-scattering surface. Even if the true pattern of anisotropies at the last-scattering surface was known, the 
finite amount of power in the CMB at small scales would still constrain lensing reconstruction. Silk damping at the 
last-scattering surface suppresses CMB power at small scales, while the finite resolution of any real experiment would 
limit the detection of any signal that is present at small scales. Lensing reconstruction fails below scales at which 
there is sufficient power, for the same reason that any remapping is indistinguishable given a uniform background. 

Here, we consider a third source of noise neglected in previous studies. The filtered temperature-squared map is an 
unbiased estimator for the lensing potential in the approximation that a correlation between two given temperature 
modes is induced only by the single lensing mode whose wavevector is the sum of that of the two temperature modes. 
In actuality, any combination of two or more lensing modes whose wavevectors sum to this total induce correlations 
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between the two temperature modes. There are many such combinations, but since these correlations add incoherently 
we do not expect a systematic bias. Nonetheless, for estimators of each individual lensing mode we must use our 
knowledge of other lensing modes to subtract off this unwanted bias. This is an iterative process, and since our 
knowledge of the lensing map is imperfect it induces noise in lensing reconstruction. We calculate this additional 
variance for various estimators constructed from CMB temperature and polarization maps, and show how it compares 
to the dominant noise sources for an all-sky CMB experiment with a noise-equivalent temperature of 1 /iK y/sec. 
Since the lcnsing-potential power spectrum is a measure of the theoretical uncertainty with which we can predict the 
value of a given lensing mode, this noise associated with lensing reconstruction causes a systematic overestimation of 
the lensing-potential power spectrum. This systematic bias must be accounted for in order to compare observations 
with theoretical predictions. 

This paper is organized as follows. In § II, we define the formalism we will use to explore the effects of weak lensing 
on the CMB. The Taylor expansion of the lensed CMB map in gradients of the lensing potential is given in both real 
and Fourier space, and the power spectra and trispectra of various components of the CMB temperature map are 
listed for later use. In § III, we show that the Fourier modes of the temperature-squared map when properly filtered 
can serve as estimators for the Fourier modes of the displacement map with the same wave vector. Using the power 
spectrum and trispectrum given in the preceding Section, we calculate the variance associated with this estimator 
including a new component neglected in previous studies. This variance is evaluated numerically using the currently 
favored ACDM cosmological model with baryon density Ot, = 0.05, matter density Q m = 0.35, cosmological constant 
density JIa = 0.65, Hubble parameter h~0.65, and power-spectrum amplitude erg = 0.9. We then use the displacement 
estimator for individual Fourier modes to construct an unbiased estimator for the lensing-potential power spectrum 
in § IV, and calculate the variance and covariance associated with this estimator. A few concluding remarks about 
the implications of our work for future studies are given in § V. The Appendix contains useful formulae related to 
additional estimators of lensing based on polarization and a combination of temperature and polarization. 



II. WEAK LENSING OF THE CMB 

We consider weak lensing under the flat-sky approximation following Refs. [7,8]. As discussed before [8,9], weak 
lensing deflects the path of CMB photons resulting in a remapping of the observed temperature pattern on the sky, 

0(n) = e[n + V0(n)] 

« 6(n) + V^(n)V J e(fi) + iv^(n)V^(n)V i V J '6(n) + . . . (1) 

where 8(n) is the unlensed primary component of the CMB in a direction n at the last scattering surface. The 
observed, gravitationally-lensed temperature map 8(n) in direction n is that of the unlensed map in direction n + 
V(/>(n) where V0(n) represents the lensing deflection angle or displacement map. Although a real CMB map will 
include secondary contributions such as the SZ effect [10], we assume that such effects can be distinguished by their 
frequency dependence [11]. They will not be further considered in this paper. A noise component denoted by 8 n (n) 
due to finite experimental sensitivity must be included as well. Thus the total observed CMB anisotropy will be 

e^n) = e(n) + e n (n). 

Taking the Fourier transform of the lensed map 8(n) under the flat-sky approximation, we write 

6(1) = J ie(n)e- M 

= ©(i)-/ ^e(i')i(i,i'), (2) 

where 

L(l, I') ee 0(1 - I') [(1 - 1') ■ 1'] + I f |^<X1") (3) 

X(j){l - 1' - 1") (1" • 1') [(1" + 1' - 1) • 1'] + . . . . 

CMB correlations in Fourier space can be described in terms of a power spectrum and trispectrum as defined in 
the usual manner, 
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(e*(li)e 4 (l 2 )) = (27r) 2 <5D(li+l 2 )C?, 

(e'(ii) . . . e i (u)) c = (2^) 2 «5 D (i 1 + 1 2 + 1 3 + uWh, h, i 3 , U) , 

(4) 

where the angle brackets denote ensemble averages over possible realizations of the primordial CMB, large-scale 
structure between the observer and the surface of last-scattering, and instrumental noise. The subscript c denotes 
the connected part of the four-point function and the superscript i denotes the temperature map being considered 
(6*, 9, or 9 n ). The lensing-potential power spectrum can be defined analogously, 

(0(l)0(l'))LSS-(27r) 2 fe(l + l')C^, (5) 

where here the angle brackets denote an average over all realizations of the large-scale structure. We make the 
assumption that primordial fluctuations at the last-scattering surface are Gaussian. Gaussian statistics are fully 
described by a power spectrum; the Gaussian four-point correlator, (0(li) . . . 9(Li)) c is zero. The instrumental noise 
O n is also assumed to be Gaussian, as is the lensing potential <f>. This second assumption is justified because the 
dominant contributions to the lensing potential come from intermediate redshifts 1 < z < 3 at which linear theory 
holds. 

Using these definitions, we can calculate the anticipated power spectrum and trispectrum of the observed CMB 
map. Because the instrumental noise is uncorrelated with the signal, the power spectrum of the observed map is the 
sum of signal and noise power spectra, 

Cf et - Cf 8 + Cf 0n . (6) 
The power spectrum of the noise component is given by: 

rieOn f „„-lJ 2 cr? (y\ 
= JskyW e » , (7) 

where / s k y is the fraction of the sky surveyed, w^ 1 is the variance per unit area on the sky, and er& = #/V81n2 is 
the effective beamwidth of the instrument expressed in terms of its full-width half-maximum resolution 9. A CMB 
experiment that spends a time t p i x examining each of -/V p ; x pixels with detectors of sensitivity s will have a variance 
per unit area w" 1 — 47r(s/TcMB) 2 /(ipixA r p i x ) [12]. The power spectrum of the lensed CMB can be determined by 
inserting Eq. (2) into Eq. (4) as discussed in [8], 
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cf e + 1 ^c^xpr^-^-hf. (8) 



This result is given to linear order in the lensing-potential power spectrum Cf . Lensing neither creates nor destroys 
power in the CMB, but merely shifts the scales on which it occurs as seen by the fact that 



-2_f d 2 \ - ee _ r d 2 \ ee _ 2 

-J W? 1 -J W 1 " 



(9) 



The observed CMB trispectrum can be calculated in a similar manner; under our assumptions of Gaussian instrumental 
noise and no secondary anisotropics the trispectrum of the lensed component 9 is the sole contribution to the total 
observed trispectrum, 



7^(11,12,13,14) =f e (l 1 ,l 2 ,l 3 ,l 4 ) 

Cg +hl [(h + 1 3 ) • l 3 ][(li + 1 3 ) • 1 4 ] + Cf* +hl [(h + 1 3 ) • ls][(l 2 + 1 3 ) • U] 



+ Perm. . 

(10) 



The term shown above is manifestly symmetric under the interchange li <-> I2, while the "+ Perm." represents five 
additional terms identical in form but with the replacement of (li, I2) and (13,14) with the other five combinations 
of pairs. The total trispectrum is symmetric under the interchange of any given pair as one would expect. Having 
established a formalism within which to analyze weak lensing, we now consider the problem of reconstructing the 
lensing potential from an observed CMB temperature map. 
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III. LENSING-POTENTIAL ESTIMATORS 



In this Section, we examine lensing reconstruction following the approach of Ref. [6] , largely adopting their notation 
as well. The only important difference in notation is that we use to denote the lensed temperature field and for 
the unlenscd field following Ref. [8] and most recent papers. Ref. [6] uses the opposite convention. For 1 ^ — 1' and to 
linear order in </>, 

(0 t (l)0 t (l / ))cMB = /ee(l,l')^(L), (") 

where 

/ee(l, I') = Cf e (L • 1) + Gf e (L ■ I') , (12) 

and L = 1 + 1'. Note that ( )cmb differs from the unmarked ( ) that first appeared in Eq. (4) in that it denotes 
an ensemble average only over different Gaussian realizations of the primordial CMB and instrument noise; a fixed 
realization of the large-scale structure is assumed. For the purposes of estimating the large-scale structure actually 
realized in our observable universe, this is the appropriate average to take to ensure that our estimators are truly 
unbiased for a typical realization of the primordial CMB. When calculating the noise associated with lensing-potential 
estimators and again for power spectrum estimation in § IV, we will return to the full unmarked ensemble average. 
Eq. (11), an immediate consqucncc of Eq. (2), suggests that a temperature-squared map appropriately filtered in 
Fourier space can serve as an estimator for the deflection field d(L) = zL0(L). Hu and Okamoto define five differ- 
ent estimators for the deflection field constructed from various combinations of the temperature and polarization; 
we discuss the temperature-squared estimator in this Section and relegate the analogous formulae for polarization 
estimators to the Appendix. The minimum-variance temperature-squared estimator derived in Ref. [6] is 

dee(L) ee iLA ®f^ J ^_@t (li)e t (l2)jFee(li)l2)) (13) 

where 

F e e(h,h) ee ^eeigeet ■ ( 14 ) 



Aee(L) = L 2 



' Hl r/ee(li, 12)^88(11,12) 



(15) 



(27T)*- 

and 1 2 = L li. Substitution of Eqs. (11), (12), (14), and (15) into Eq. (13) shows the desired result, 

(d ee (L))cMB =d(L), (16) 

namely that dee(L) is indeed an unbiased estimator for the deflection field in Fourier space. We now proceed to 
calculate the variance of this estimator. At first, we assume a complete knowledge of all lensing modes not examined 
by this estimator. In that case, we find 

(d@@(L) • d e(L')) C MB - (dQ (L)) CM B ■ (dee(L'))cMB = (L • L') e ®J - e ® /2 — - 

x J J |^(0 t (-li)0 t (-l 2 )0 t (li')e t (l 2 '))cMBFee(li,l2)Fee(li',l2')-d*(L).d(L'), 

(17) 

where I2' = L' — li'. Evaluating the four-point function in the integrand of Eq. (17) to second order in the lensing 
field, we obtain 
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FIG. 1. The two quadrilaterals consistent with the constraint L — L' = h + 12 — li' — h' = for the variance of the estimator 
dee(L)- The lensing modes (j>(L), 4>(h — li'), and 0(li — h'), depicted as diagonals in the above quadrilaterals, induce 
non-Gaussian couplings between the modes of the observed temperature map represented as sides of these quadrilaterals. They 
lead to the three groups of linear terms appearing in Eq. (18). 



(e t (-i 1 )e t (-i 2 )e t (i 1 ')e t (i 2 '))cMB 



(c ee + c ee n) (27r) 2 fo(L) + 0(_L)/ ee (li, 1 2 ) - / ^2 C v @ <f>(-h l')<K~h + 1') [l' ■ (h + 1')] 

t 1 ' • (l2 ~ ~ \ I -0ji^w- L - nj^rai • n [h ■ & + o] + • n [i 2 • (l + n] } 

(Cf e + Qf n ) (2^) 2 ,5 D (L') + 0(L')/ee(li', b') - / ^^■C^ e V(li' - l'Mfe' + 1') [l' ■ (U' - 1')] 
t 1 ' ' (l2 ' + + \ j ^ HmV V) { C % @{h ' ' 10 1 11 ' ' (1 ' ~ L ' } ] + 8(l2 ' ' 10 t l2 ' ' (1 ' ~ L ' } ] I 



+Pcrm. 
(18) 



The terms given explicitly in Eq. (18) correspond to the correlations between ©*(— li) and ©*(— 12) and those between 
t (li') and Q t (h')- The "+ Perm." stands for two additional terms, identical in form, arising from the pairings 
(0 t (-l 1 )0 t (l 1 ')>CMB(0 t (-l2)0 t (l2 , )>CMB and (0 t (-li)0 t (l2 , )>CMB(0 t (-l 2 )0 t (li'))cMB. This expression indicates 
how uncertainty in the CMB at the last-scattering surface propagates into uncertainty in lensing reconstruction for a 
particular realization </>(L) of the large-scale structure. To linear order in Eq. (18), correlations between the modes 
©*(— li), ©*(— 12), t (li'), and ©'(l/) are induced by those lensing modes whose wavevectors are the sums of any pair 
of wavevectors of these four modes. These lensing modes are precisely those forming the diagonals of the quadrilaterals 
depicted in Fig. 1. In practice, we do not know the large-scale structure between us and the last-scattering surface, 
so we assume a variance given by Eq. (5) with a model-dependent power spectrum Cff. We must average Eq. (17) 
over different realizations of the large-scale structure (denoted by ( )lss) to obtain the total expected variance of 
our estimator, 

( d ee( L ) ' d e e(L'))cMB - (d<L )e (L)) CM B ■ (dee(L'))cMB) Lgs = (L • \,') Ae ®^> - 

X S Wf! 7^( et (- 1 i) 0t (-l2)0 t (li')0 t (l2')>^ee(li,l 2 )^ee(li',l 2 ')-(2^) 2 <5D(L-L')Cf , 



where 



(d*(L) • d(L')) = (2tt) 2 <5 d (L - I/)CT = (2tt) 2 ,5 d (L - L')L 2 C£ 



(19) 



(20) 



The assumption that the lensing potential is Gaussian imposes the constraint L — L' = li + 1 2 — li — I2 = which 
closes the quadrilaterals of Fig. 1. The average of the four-point correlation function in Eq. (19) can be calculated 
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by further averaging Eq. (18) over the large-scale structure. Terms linear in the lensing field vanish when averaged 
over different realizations of the large-scale structure. Quadratic terms in the lensing field arise as products either of 
two linear terms or of a zeroth and second-order term. Averages over the product of two linear terms produce the 
connected part of the four-point correlation function, the trispectrum defined in Eq. (4). Averages over the product 
of a zeroth and second-order term have no connected portion, but instead furnish an implicit dependence on Cf/ 
in the total observed power spectrum of Eq. (6). The final result of averaging over the large-scale structure can be 
expressed in terms of the observed power spectrum and trispectrum, 

(6 t (-l 1 )6 t (-l 2 )e t (l 1 ')e t (l 2 ')> - (2tt) 4 [cf^Cf e %(L)ML') + c/f^Qf^Mii' - h)Su(h' - h) 

+5 u (l 2 ' - li)Mli' -h)} + (27r)- 2 T t (-l 1 , -1 2 , 1/, 1 2 >d(L - L')] , 

(21) 

where the trispectrum can written in terms of /ee(li>l2) & s 

T\h, h, h, U) = Cjf +la| / ee (li, l 2 )/ee(l3, U) + qt?+i 3 |/ee(li, h)fee(h, h) + Cff+u | fee (h, h)fee(h, h) ■ (22) 

This form of the trispectrum is consistent with that of Eq. (10) given directly in terms of power spectra. Since L^O, 
c>d(L) = and the first term of Eq. (21) vanishes. The remaining two terms containing pairs of delta functions, 
inserted into Eq. (19), yield the dominant contribution to the variance, 

(d* @e (L) ■ dee(L')) C MB - (d^ e (L)) CM B ■ (d e e(L'))cM B ) - (2^) 2 <5 D (L - L')[N^ ee (L) + ...], (23) 

/ LSS 

where Nqq ee (£) = Aq®(L). Notice that N Q ° @ @& (L) is zeroth order in the lensing potential (ft; it depends on the 
lensing potential power spectrum Cf^ only implicitly though the total observed power spectrum C ; eet . The ellipsis 
represents terms of higher order in Cf that we now proceed to calculate. These terms arise from the trispectrum 
tcrmofEq. (21) after ((d0@(L))cMB ■ (dee(L'))cMB)LSS = (27t) 2 <5d(L — V)Cff is removed. Substituting these results 
into Eq. (19), we find that to first order in Cf^ , 

(d^ e (L) • dee(L'))cMB - (d£ e (L)) C MB ■ (d@@ (L'))cmb\ qq - (2^) 2 <5 D (L - L')[JV(°> (L) + N e % @@ (L)] , (24) 

where N^l <ee (L) is given by, 

(i) A% B (L) f d 2 h f d% 

N ee,eey L ) - — J2 — J J^yl J j^T2 F ee{h,h)F @ e{h ,h ) 

x {Cf^feei-h, h')fee(-h, W) + C++_ W \fee(-li, h')fee(-h, h')} • 

(25) 

The first-order contribution to the noise Nqq @& (L) involves integrals over the lensing-potential power spectrum, 
and thus probes lensing modes with wave vectors different from that of the estimator dee(L). It can be interpreted 
physically as interference from these other modes in the determination of the mode d(L) being estimated. The filter 

Fqq(1i,1 2 ) was chosen to optimize the signal-to-noise ratio in the absence of the first-order contribution Nqq ee (L); 

it is no longer an optimal filter once this additional noise is taken into account. As long as Nq @ && (L) <§C N @ ° @ ee (L), 
the noise reduction that can be attained by re-optimizing our filter will not be significant. Formulas analagous to 
those presented here relevant to the construction of estimators using polarization data are given in the Appendix. 

The significance of Nqq @@(L) for two different experiments is shown in Fig. 2 using the currently favored ACDM 
cosmological model with baryon density = 0.05, matter density fl m = 0.35, cosmological constant density 17a = 
0.65, Hubble parameter h~0.65, and power-spectrum amplitude er 8 = 0.9. The Planck experiment is equivalent to a 
one-year, full-sky survey with temperature and polarization sensitivities of 12.42 and 26.02 fiKy/sec respectively and 
resolution = 7.0 arcminutes as described in § II. The reference experiment has the same resolution but superior 
sensitivities of 0.46 and 0.65 /iK^/sec for temperature and polarization. These estimates of experimental parameters 
are identical to those given for the Planck and reference experiments of Ref. [6] . The QE and EE estimators have noise 
power spectra intermediate to those of the 96 and EB estimators, while the QB estimator has substantially higher 
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FIG. 2. Variances with which individual modes d(L) of the deflection field can be reconstructed by the Planck and reference 
experiments described in the text. The solid curves are the power spectra Cf d anticipated for our ACDM cosmological model. 
The upper and lower dashed curves are the zeroth and first-order noise power spectra Nqq ee {L) and Nqq e@ {L) respectively 
for the temperature-based estimator dee(L), while the dotted curves are the corresponding noise variances for d_Es(L). A 
mode d(L) cannot be reconstructed with signal-to- noise greater than unity when C[ d < ee (L) + @0 (L). 



noise because the primordial CMB lacks true B-modes in the absence of inflationary gravitational waves. We see that 
for Planck, with its comparatively inferior polarization sensitivity, the 96 estimator will be best although it will be 
unable to detect individual Fourier modes of the deflection field at the la level. The reference experiment, and further 
experiments with similar sensitivity and even higher resolution, should be able to push lcr detection of individual d(L) 

modes to L ~ 1000 by primarily relying on the EB estimator. In these cases the secondary noise N BB EB (L) is only 

smaller than the dominant noise N EB EB (L) by a factor of a few, whereas for higher sensitivity (noisier) experiments 
like Planck it is smaller by at least an order-of-magnitude. This illustrates an interesting point, apparent from Fig. 2, 
that the zeroth-ordcr noise Nqq &@ (L) declines dramatically with decreasing sensitivity until it becomes dominated 
by cosmic variance while the Nqq 0e (L) is largely unaffected by instrument sensitivity. The reasons for this trend are 
that instrument noise appears in Nqq @&(L) = A@®(L) through its contribution to the denominator of Fee(li, I2) as 
shown by Eqs. (14) and (15). Decreasing instrument noise raises the value of Fee(li, I2) thereby lowering Nqq ee (I/). 
By contrast instrument noise is reflected in Nqq qq{L) through its effects on both A@®(L) and i*ee(li, h) as shown 

in Eq. (25). Smaller instrument noise raises i*ee(li,l2) as before, driving Nqq ee (L) up in this case, but this is 
compensated by a decrease in A@@(L) which appears as a prefactor outside the integrals. These two effects largely 
cancel each other out, rendering Nqq qq(L) remarkably insensitive to instrument noise. 
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IV. POWER SPECTRUM ESTIMATION 



Although complete reconstruction of the deflection field d(L) can be an enormously powerful tool, such as for 
B-mode subtraction [4], for some purposes estimates of the lensing-potential power spectrum Cff are sufficient. 
This power spectrum is a model-dependent prediction of theories of large-scale structure formation, and therefore 
estimates of the power spectrum from real data could be used to test these theories as well as the consistency of other 
determinations of cosmological parameters. Furthermore, since estimates of all the modes d(L) with |L| = L can be 
combined to estimate Cff ', la detection of the power spectrum can be pushed to much higher L than can that of 
individual modes. The deflection-field estimator dee(L) derived in the preceding Section can be used to construct 
an estimator for Cff. Our first guess for an appropriate lensing-potential power spectrum estimator is 

D ^^^J aL (^ d ^-aee(-l), (26) 

where A is the area of the sky surveyed and is an annulus of radius L and width AL. We ensemble average our 
estimator over different realizations of the CMB and large-scale structure using Eq. (24) by bringing ((d@ e (L))cMB ■ 
(dee(L'))cMB)LSS from the left to the right-hand side. This yields 

{Dl) = aI?2^aZ MQ) / dH + ^ee.eeCO + <lee(0] • (27) 

The definition of the Dirac delta function, 

(2tt) 2 (5 d (1) = / dhe a A , (28) 

J A 

implies that <5d(0) = A/(2ir) 2 . Furthermore, in the limit that AL is small compared to the scales on which l 2 Cf^ + 
Nqq ee (0 + Nqq 0q (0 is varying, we can evaluate the integrand of Eq. (27) at its central value I = L and extract it 
from the integral. The integral over the annulus cancels the factor 2nLAL in the denominator, reducing Eq. (27) 
to 

(D L ) = Cf + L~ 2 [Nj» >ee {L) + Ngl @e (L)] . (29) 

Dl is indeed an estimator for the lensing-potential power spectrum Cf/, albeit a biaed one. Note that the bias in 
the power-spectrum estimator Dl is precisely the same as the variance shown in Eq. (24) with which we were able 
to determine each individual lensing mode. This is no coincidence; it reflects the fact that there are no grounds a 
priori on which to differentiate the variance with which we can reconstruct individual modes d(L) from the intrinsic 
variance Cjf of the underlying distribution from which they are drawn. To obtain an unbiased estimator to compare 
with theoretical predictions, we subtract off this unwanted reconstruction variance, 



^ee(i) + Jv4e.ee W ■ ( 30 ) 



Since Nqq qq{L) as defined in Eq. (25) itself depends on Cff, this subtraction and evaluation must be performed 

iteratively until a self-consistent solution is obtained. The variance of our estimator Cf/ can be calculated in the 
usual manner, 

= ({Cf) 2 ) (Cf) 2 = ((D L ) 2 ) (D l ) 2 . (31) 

L 

Evaluating this expression requires us to calculate 

« Di > 2 > - %(dKE?L L w? <[dee(1,,) ■ dee( - |, ' )1 [dee(l2,) ' deeH2 ' )]> • (32) 

Since d©©(L) is a quadratic estimator in the temperature map, Eq. (32) includes the following integral over the 
eight-point correlation function in Fourier space, 



8 



/[j n H , , M1 r H /i /\ j , , ,viv ^ee(^i)^ee(^) / d 2 k, f d 2 k 3 f d 2 k 5 f d 2 k 7 
([deeOi ) • dee(-li )J [d ee (l 2 ) ■ d e e(-l 2 )J) - ( ,, )2(z , )2 j J J J 

x{fee(ki,k 2 )^ee(k3,k4)^ee(k5,k 6 )Fee(k7,k 8 )(e t (k 1 )e t (k 2 )e t (k3)e t (k4)e t (k 5 )e t (k^ 

(33) 

where k 2 = l/ — ki, k 4 = — li' — k 3 , k 6 = 1 2 ' — k 5 , and k 8 = — 1 2 ' — k?. A fully general eight-point correlation function 
consists of a connected part, as well as terms proportional to the product of lower-order correlation functions. Under 
the assumption that both the primordial CMB and the lensing potential are governed by Gaussian statistics, all 
correlation functions higher than the four-point have vanishing connected parts [13]. The temperature eight-point 
correlation function will therefore be composed of three groups of terms; membership in a group being determined by 
whether the term contains zero, one, or two factors of the trispectrum. Since the trispectrum given in Eq. (10) is first 
order in the lensing-potcntial power spectrum Cjj, terms of these three groups are zeroth, first, and second order 

respectively in Cff. Combinatorics determines the number of terms in each group. There are: ^ ^2) (2) (2) (^2 
105 different ways of dividing (ki, . . . , k 8 ) into four pairs, and hence there will be 105 terms in the group containing 
no trispectra. Similar calculations reveal that there arc ^ = 210 terms in the second group and 

2l (a) (a) = ^ terms in the third group. Many terms in all three groups will vanish for the same reason that the 



v 4 / \ 4 y 

first term vanished in the four-point correlation function of Eq. (21); these terms are proportional to a Dirac delta 
function evaluated at nonzero argument. Consider now the first group of terms, those that are zeroth order in Cf^. 
The 60 nonvanishing terms in this group each contain four Dirac delta functions; they can be further segregated into 
the 12 terms that allow two of the integrals over k^ appearing in Eq. (33) to be immediately evaluated via Dirac delta 
functions, and the 48 terms that allow evaluation of three ki integrals. The first 12 terms, inserted into Eq. (33) and 
appropriately evaluated using the normalization of Eq. (15), yield 



((Dl) 2 ) 



'Ngl t@@ (L)\ - 



L 2 



2 w 1 



A 2irLAL 



+ ..., (34) 



while the remaining 48 terms give the final result to zeroth order in Ctf , 



{(Dl) 2 ) 



L 2 
+ (2tt) 



2 



(27T) 



2 



A 2ttLAL 



AL A (2ttLAL) 



f d 2 h> f d 2 h> A 2 ee (l[)A 2 ee (l' 2 ) f d 2 k, 

(35) 



where 



P(ki, k 2 , 1/, h') = fee(-h' - ki, h' - k 2 )F ee (-ki, 1 2 ' + k 1 )F ee (-k 2 , -1 2 ' + k 2 ) 
+fee(h' - ki, -1 2 ' - k 2 )F ee (-ki, -h' + k 1 )F ee (-k 2 , h' + k 2 ) 
+/ee(-l 2 ' - ki, li' - k 2 )Fee(-ki, h' + k 1 )F ee (-k 2 , -h' + k 2 ) 
+/ee(li' - ki, -1 2 ' - k 2 )F ee (-k 1 , -1/ + ki)F ee (-k 2 , h' + k 2 ) 
+/ee(li' - ki, h 1 - k 2 )F 8 e(-ki, -li' + k 1 )F ee (-k 2 , -1 2 ' + k 2 ) 
+/ ee (l 2 ' - ki.li' - k 2 )F e e(-ki,-l 2 ' + k 1 )F ee (-k 2 ,-li' + k 2 ) . 



(36) 



When {Dl) 2 is subtracted from {(Dl) 2 ) in Eq. (31), the first term in the square brackets of Eq. (35) will be eliminated. 
Minimizing the variance associated with this estimator then consists of making an optimal choice of AL(L). The first 

noise term is proportional to ^kIal ■ ^ or a survev 01 area A (2tt) 2 /A is the specific area of an individual mode in 
L space and 2-kLAL is the area in L space over which the power-spectrum estimator takes an average. This ratio 
is therefore the inverse of the number of individual dee (L) modes whose inverse variances are added to determine 



9 



the inverse variance of Cjj. It is obviously minimized by choosing (2ir) 2 /A <C 2irLAL(L). The second term, that 
involving P(ki, k 2 , li , 1 2 ), differs from the first noise term in that a Dirac delta function has been used to evaluate 
an additional integral rather than an annulus integral. Since the integrands are of the same order, we expect the 
second noise term to be suppressed relative to the first by a factor 2nLAL(L) /7rZ^ ax where l max ~ ir/9 is set by the 
resolution 9 of the survey. Under the conservative assumption LAL(L) <C 1/9 2 , namely that we are probing scales 
well above our resolution, this term is assured to be small. We neglect such terms for the remainder of this paper. If 
we insert the portions of the eight-point correlation function that are first and second order in Cff into Eq. (33) and 
evaluate using Eq. (25), we find 



{(Dl) 2 ) = (cf + N$l ee (L) + Ngl e@ (L)) ' 



1 + 2 



(2^) 2 



1 



A 2nLAL 



(37) 



and 



2 

O g<f,4> 



= 2 



(2tt) 2 



1 



A 2ttLAL 



(38) 



This result agrees with that given in Ref. [6] after subtracting our newly derived term Nqq @@(L). 



The term Nqq @q(L) and corresponding terms for polarization-based estimators have two principal effects on 
power spectrum estimation. As shown in Eq. (38), they provide a fractional contribution to the variance of roughly 
(L)/Cff when Cjf dominates the variance as in the reference experiment in the right-hand panel of Fig. 2. 



(i) 

ee.ee 



2N ( 

For the ACDM cosmological model considered here this represents an increase of 5 — 15% in the variance of the EB 
estimator for L < 1000. More importantly, Nqq q@(L) acts as a bias for the naive estimator Dl as shown by Eq. (29). 

If this bias is not calculated and subtracted iteratively to form the unbiased estimator Cf/ as in Eq. (30), Cff will 
be systematically overestimated by 5 — 10% at low L and by increasingly larger amounts at L > 100 as the signal 

L(L + l)C"[ d /2ir begins to plummet while Nqq q@{L) remains comparatively flat. 

Having evaluated the variance of our estimator Cff, we consider whether this estimator has a substantial covariance 



a et , e , f ee <(cf - Ct«)(Cit - C+*)) = (D L D L ,) (D L )(D L .) . 
The estimator Dl as defined in Eq. (26) implies that: 



(39) 



{D L D L >) 



2tt 



{'ALL' J LALL' 



17au J aL Wf I , Wf <[dee(ll,) ' dee( - ll/} ] [ dee(l2 ' } ' d0e (- l2 ')]> ' ( 4 °) 



which can be evaluated using the same integral over the eight-point correlation function described in Eq. (33). 
Whereas 60 of the 105 zeroth-order terms in Cf/ coming from this equation were nonvanishing for the variance, for 
the covariance only 52 terms are nonzero provided that the widths AL and AL' are chosen so that the annuli o_l and 
a,L' do not overlap. This leads to a result analogous to Eq. (35), 



{D L D L >) = - 7 



(o) 



Li 1 



L' 2 



+ 

2l / 



AL 3 ALL 



'3AL' J aL (2tt) 2 L lI (27T 



d 2 h' Al e (l[)A 2 ee (l> 2 ) f d 2 k, 



2 (i'i) 2 ^) 2 



f d 2 k 

J (2^/ ee ( kl ' k2 )' P ( kl ' k2 ' 11 '' 12 ') ' 



(41) 



Note that the eight terms missing from the covariance when compared to the variance have altered the first term of 
Eq. (41), and that it was precisely these terms that provided the dominant contribution to the variance of Eq. (38) 
when AL(L) was chosen appropriately. We therefore find that to zcroth order in Cf/, the covariance is given by 

2 



«J ee (I)C,9e(i') 



(0) 



d 2 h' 



ALALL'AL' 



L 2 



L' 2 



d 2 h' 



aL (27T) 2 J ar , (27T) 2 J (27T) 2 



d 2 k! 



/ee(ki,k 2 )P(k 1 ,k 2 ,l 1 ',l 2 ' 



(42) 
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FIG. 3. The ratio Rll' for Planck as a function of L for fixed values of V . The solid curves correspond to L' = 3, 30, 300 as- 
cending from bottom to top, while the long-dashed and short-dashed curves correspond to L' = 7, 70, 700 and L' = 10, 100, 1000 
respectively, again with curves in each sequence appearing from bottom to top in the figure. 



where we have extracted the Aee(^) from the annular integrals since they are slowly varying over the widths AL 
and AL'. For the same reasons that terms of this form were a subdominant contribution to the variance as discussed 
previously, we expect the covariance to be suppressed as well. If we define the ratio 



R 



LL> 



(43) 



q4><P 



we can quantify this suppression. The triple integral of Eq. (42), appearing in the numerator of Rll 1 , involves 
integration over annuli with radii L and L' and one integration over all Fourier space. The triple integrals in the 
variances d 1 -^ and cr-^ appearing in the denominator of Rll 1 each consist of a single integration over an annulus 

of radius L and L' respectively and two integrations over all Fourier space. If we make the crude assumption that the 
integrand is constant, the ratio Rll' will simply be the ratio of these areas, 



Rll 1 — ^{2itLAL){2itL'AL')/itI 



2 

max' 



(44) 



The ratio Rll 1 is evaluated numerically for Planck in Fig. 3 as a function of L for various fixed values of L' . The 
estimators Cf^ and Cff were chosen such that AL = AL' = 1, while integrals over Fourier space were cut off at 
( ma![ = 5000. Substituting these values into Eq. (44), we expect Rll 1 — 8.0 x 10~ 8 V LL' . This crude estimate is 
surprisingly close to the numerically obtained results of Fig. 3; in particular the slope of the curves is approximately 1 /2 
on this log-log plot. Even at L ~ 1000, Rll 1 ^ 4.0 x 10~ 4 suggesting that covariance in power-spectrum estimation 
can safely be neglected for Planck. The estimate of Eq. (44) implicitly depends on the experimental resolution 9 
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because the integrand appearing in expressions for the variance and covariance decreases rapidly for L > Z max ~ tt/6. 
For future experiments with better resolution than Planck, ^ max will be higher implying by Eq. (44) that covariance 
will be even more negligible. 

V. DISCUSSION 

Weak gravitational lensing induces non-Gaussian correlations between modes of the observed CMB temperature 
map as shown in Eq. (11). These correlations, and assumptions about the Gaussian nature of the primordial CMB, can 
be used to construct several temperature and polarization-based estimators of the Fourier modes d(L) of the deflection 
field. This procedure was outlined in Ref. [6] , however in calculating the noise associated with this reconstruction, an 
assumption was made that the observed temperature map was Gaussian. In the presence of lensing this assumption 
is invalid; when calculating the variance of quadratic estimators all permutations of the observed trispectrum must 
be taken into account. One such permutation reflects the desired correlation making our estimator sensitive to d(L), 
but the remaining two permutations induce additional variance proportional to the lensing-potcntial power spectrum 
Cf/. While subdominant, this variance will become increasing significant for future experiments as shown in the 

right-hand panel of Fig. 2. Since the power spectrum Cf^ is itself a measure of uncertainty in the deflection field, 
this additional variance in lensing reconstruction acts as a bias during power-spectrum estimation because there is 
no a priori way to distinguish it from the intrinsic variance of the underlying distribution. Our calculation of the 
dependence of this variance on Cf/ allows it in principle to be subtracted iteratively, which will prevent a systematic 
5 — 10% overestimate of Cff at low L. 

We close by considering several possible observational obstacles to the scheme for lensing reconstruction and power- 
spectrum estimation presented above. One hindrance is secondary contributions to the CMB such as the SZ and 
ISW effects. These effects increase the total temperature power spectrum appearing in the denominator of the 
optimum filter Fee(li> I2) of Eq. (14) as would additional instrumental noise. They also correlate with the large-scale 
structure at low redshifts inducing further non-Gaussian couplings and additional variance to lensing reconstruction. 
Fortunately for our purposes the frequency dependence of the thermal SZ effect differs from that of a blackbody. It can 
therefore be separated in principle from the lensed primordial CMB by an experiment with several frequency channels 
[11]. The ISW effect cannot be removed in this manner, but is too small to significantly inhibit lensing reconstruction. 
Polarization-dependent secondary effects are expected to appear at higher orders in the density contrast [14], and 
we therefore anticipate that they will not make a contribution at the levels considered here. A potentially more 
serious problem is that of galactic foregrounds, which though uncorrelated with the lensing signal may be substantial 
at certain frequencies. Significant polarization has also been observed in some of these sources [15]. We hope to 
understand and minimize the effects of galactic foregrounds in future work, and to pursue further refinements of 
lensing reconstruction. 

ACKNOWLEDGMENTS 

We thank Wayne Hu for useful discussions. This work was supported in part by NASA NAG5-11985 and DoE DE- 
FG03-92-ER40701. Kesden acknowledges the support of an NSF Graduate Fellowship and AC acknowledges support 
from the Sherman Fairchild Foundation. 



APPENDIX A: POLARIZATION-BASED ESTIMATORS 

Here we provide the appropriate formulas for deriving the variance associated with polarization-based estimators 
of the deflection field d(L). The CMB polarization can be decomposed into E and B-modes [16]. These modes are 
mixed by weak lensing such that to linear order in <j>{L), 

E{\) = E{\) - J ^ [E(h) cos 2^ - Vl ) - B(h) sin 2^ - pi)] 0(1 - Ii) [(1 - li) • li] , 

B(l) = B(l) - J [E (h) sin 2( Vll - W ) + B(h) cos 2( Wl - &)] 0(1 - h) [(1 - li) • li] . (Al) 

We can exploit the sensitivity of the polarization modes to the lensing potential to construct lensing estimators from 
quadratic combinations of polarization modes. We generalize Eq. (11) to arbitrary combinations {X, X'} of O, E, 
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and B-modes as first derived in Ref. [6], 

{X* (l)X'* (1'))cmb - I xx' (1, , (A2) 

where 

f @E (l, 1') = Cf E cos2(^ - W )(L • 1) + Cf E (L ■ 1') = / B0 (l', 1) , 

/e B (l, 1') - Cf E sin 2fa - W )(L • 1) = / B e(l', 1) , 

1') = [Cf £ (L • 1) + C E £ (L • I')] cos2(^, - W ) , 

1') = [C EE (L ■ 1) + C BB (L • I')] sin2(^ - ^) = / B£ (I', 1) , 

f BB (l, 1') = [Q BB (L • 1) + C BB (L • I')] cos2(^ - W ) . 

(A3) 

In deriving these results we used parity considerations to demand Cf B = C EB = 0. Using these relations, we follow 
the approach of Ref. [6] to derive symmetric lensing estimators 



iLA xx (L) f d% 

and 



d.w(L) z ' f ^X I \ I / , , I .1 . (A.4) 



iLA xx ,{L) f d\ 1 

dxx'[L>) = 



L 2 



I Wf\ [^(li)*^) + X'\h)X\h)] F xx ,{\ u h) . (A5) 



We have explicitly symmetrized our estimators for X ^ X' to simplify the form of the optimal filter. The normalization 
bias of the estimators is removed by choosing 



A XX (L) = L 

and 

A XX ,(L) = L 2 



/cPli 
j^y}fxx{h,h)Fxx{h,h) 



J J0?\ [fxx ' {luh) + fxx ' (l2 ' Fxx ' {h ' h) 



(A6) 



(A7) 



The minimum- variance filters Fxx'(h,h) for the various cases {X, X'} are given by 

p /lis /eg(li, h) + feE(h,h) 

r@E(n, 12) — ^eet^BBt 1 2C eEt C 0Et + c EEt C e0t ' 



F &B (h,h) 
Fee (h,h) 
Feb (h,h) 
Fbb (h,h) 



Q@8t(jEEt 1 (jEEt(j&et 

f EE (h,h) 
2(jEEt(jEEt ' 

fEB(h,h) + f EB (h,h) 
rjEEtfjBBt 1 (jBBtfjEEt 

f BB (h,h) 
2(jBBtr<BBt ' 



(A8) 



Using these optimal filters for the estimators defined in Eqs. (A4) and (A5), we can calculate the variances for these 
estimators in a fashion entirely analogoug to Eq. (19), 

( d xx( L ) ' d X x(L'))cMB - (dx X (L)) CM B • (dxx(L'))cMB^ Lgs = XX ^ - XX j — - 
X /(S^/ |^(^ t (- 1 l)A t (-l2)A t (l/)X t (l 2 0)fe(ll,l2)^ X (ll^l20-(2^) 2 5D(L-L')Cf , 

(A9) 
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'/H* (T\ H ri'U /j* /j , T ,« \ ^x'W A XX ,(L') f d% f d% 
^ d xx>\ L ) ■ axx'(L, )) cmb - {dxx'{ L ))CMB ■ [axx'K^ ))cmb / > Lgg = ^ ^ y T^p- y 



(2tt) 2 y (2tt) 2 

xi([X t (-l 1 )X' t (-l 2 )+^' t (-li)X t (-l 2 )] [^(l/^^O+X'^hOX*^')])^^!!,^)^^'.^') 
-(27r) 2 <5 D (L-L')Cf • 

(AlO) 

As for that of the temperature estimator, these variance will consist of zeroth-order terms in Cff, N xx XX (L) = 
Axx(L) and N xx , XX ,(L) = Axx'(L), and first-order terms, 



N'xkxxiL) = / / —^F X x(h,h)Fxx(h\h') 



(i) fT ,_A xx (L) f dH, f d% 

(27T) 2 J (2^) 2 

{c' | tf_ li , | /x X (~li,liO/xx(-l 2 ,l20 + C' | tf-i 2 '|/^(- 1 i' l2 ')/^(- l2 J 1 i')} 



(All) 



V (D m A\ X ,{L) f d\ f d% 

N XX>,XX>(L) = ^ J (2^)2 y ^2fc'(ll,l2)ta'(ll ,1 2 ) 

xJ{qtf-l 1 n[/^(- 1 l' 1 l')/x'X'(-l2,l2') + /xX'(-ll,ll')/^x(-l 2 ,l2') 
+/X'x(-ll, ll')/xX'(-l2, 1 2 ') + /jf'JC'(-ll, ll')/xjc(-l2, 1 2 ')] 
+qtf-l 2 '|[/^(-ll,l2 / )/x'x(-l2,ll') + /xx(-ll,l2')/x'X'(-l2,ll') 

+f X >X> (-h,h')fxx(-h,h') + fx> X (-h,h')fxX' (-12, h')] } . 



(A12) 



The filters given in Eq. (A8) are no longer optimal in the presence of this additional noise, but the difference between 
these filters and the optimal filters should be negligible provided that N xx XX {L) <C N xx XX (L), N xx , XX ,(L) <C 

N xx , XX ,(L). For the purposes of power-spectrum estimation, the terms N xx XX (L) and N xx , XX ,(L) arc not 
only an additional contribution to the variance, but are also a systematic bias if not subtracted iteratively following 
Eq. (29). 

A final point to consider is that the six different estimators d XX i(L) defined in this paper are not independent, as 
they are constructed from only three distinct maps. The covariance matrix for the six estimators will therefore not 
be diagonal, and this needs to be taken into account if the estimators are to be linearly combined to produce a single 
minimum- variance estimator. The off-diagonal elements of the covariance matrix can be evaluated in a straightforward 
manner involving pairs of double integrals similar to those of Eqs. (A9) and (AlO). 
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